Electroweak baryogenesis in the standard model with strong hypermagnetic fields 
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We show that in the presence of large scale primordial hypermagnetic fields, it is possible to generate 
a large amount of CP violation to explain the baryon to entropy ratio during the electroweak phase 
transition within the standard model. The mechanism responsible for the existence of a CP violating 
asymmetry is the chiral nature of the fermion coupling to the background field in the symmetric 
phase which can be used to construct two-fermion interference processes in analogy to the Bohm- 
Aharanov effect. We estimate that for strong hypermagnetic fields By = (0.3 — 0.5) T 2 the baryon 
to entropy ratio can be ps/s = (3 — 6) x 10 -11 for slowly expanding bubble walls. 
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The understanding of the mechanism that generates 
the observed excess of baryons over antibaryons in the 
universe represents one of the most challenging problems 
for particle physics as applied to cosmology. Any theory 
that aims to explain such excess has to meet the three 
well know Sakharov conditions Q, namely: (1) Existence 
of interactions that violate baryon number; (2) C and 
CP violation; (3) departure from thermal equilibrium. 
It is well known that the above conditions are met in 
the standard model (SM) of electroweak interactions if 
the electroweak phase transition (EWPT) is of first or- 
der. This raised the interesting possibility that the cos- 
mological phase transition responsible for particle mass 
generation -that took place at temperatures of order 100 
GeV- could also explain the generation of baryon num- 
ber. Consequently, a great deal of effort was devoted to 
explore such a possibility Q. Nowadays, the consensus 
however is that the minimal SM as such, cannot explain 
the observed baryon number. The reason is that the 
EWPT turns out to be only too weakly first order which 
in turn implies that any baryon asymmetry generated at 
the phase transition was erased by the same mechanism 
that produced it, i.e. sphaleron induced transitions ||. 
Moreover, the amount of CP violation coming from the 
CKM matrix alone cannot account by itself for the ob- 
served asymmetry, given that its effect shows up in the 
coupling of the Higgs with fermions at a high perturba- 
tive order , giving rise to a baryon to entropy ratio at 
least ten orders of magnitude smaller than the observed 
one. 

Nevertheless, it has been recently pointed out that, 
provided enough CP violation exists, the above scenario 
could significantly change in the presence of large scale 
primordial magnetic fields |||| (see however Ref. 0), 
which can be responsible for a stronger first order EWPT. 
The situation is similar to a type I superconductor where 
the presence of an external magnetic field modifies the 
nature of the superconducting phase transition due to 
the Meissncr effect. 



Though the origin and nature of these primordial fields 
is a subject of current research, its existence prior to the 
EWPT epoch can certainly not be ruled out j^] . The only 
constraint on their strength comes from considerations 
on the measured anisotropies of the cosmic microwave 
background and on nucleosynthesis calculations 

Recall that for temperatures above the EWPT, the 
SU(2)xU(l)y symmetry is restored and the propagat- 
ing, non-screened vector modes that represent a magnetic 
field correspond to the U(l)y hypercharge group. Thus, 
in the unbroken phase, any primordial magnetic fields be- 
long to the hypercharge group instead of to the U(l) em 
group and are therefore properly called hypermagnetic 
fields. 

In this paper we show that the existence of such pri- 
mordial hypermagnetic fields also provides a mechanism 
to produce a large enough amount of CP violation dur- 
ing the EWPT to possibly explain the observed baryon 
to entropy ratio in the SM. This can happen during the 
reflection of fermions off the true vacuum bubbles nu- 
cleated during the phase transition through an interfer- 
ence process equivalent to the Bohm-Aharanov effect, 
given that in the unbroken phase, fermions couple chi- 
rally to hypermagnetic fields with the hypercharge. The 
chiral nature of this coupling implies that it is possi- 
ble to build a CP violating asymmetry dissociated from 
non-conserving baryon number processes that can then 
be converted to baryon number in the unbroken phase 
where sphaleron induced transitions are taking place with 
a large rate. The existence of such asymmetry provides 
a bias for baryon over antibaryon production. In the ab- 
sence of primordial magnetic fields, this mechanism has 
been proposed and studied in Refs. fl(i|-|l"i| in extensions 
of the SM. 

To describe the EWPT, we start by writing the effec- 
tive, finite temperature Higgs potential which, including 
all the one-loop effects and ring diagrams, looks like 



^eff (0,T) = |(T 2 ~T c 2 ) 
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where <\> = v / 2( < i )tc f ) ) 1 / 2 is the strength of the SU(2) Higgs 
doublet $ whose vacuum expectation value v is given by 



(2) 



The parameters 7, 5 and A have been computed pertur- 
batively to one loop and can be expressed in terms of v, 
the SU(2) gauge boson masses and the top mass. Their 
explicit expressions can be found elsewhere (see for exam- 
ple Ref. ||). 5 is the parameter responsible for the first 
order nature of the phase transition. It is the parame- 
ter that gets enhanced in the presence of hypermagnetic 
fields. T c is the critical temperature at which spinodal 
decomposition proceeds. 

We can write the effective potential in a more trans- 
parent form fl3| by introducing the dimensionless tem- 
perature 1? and the dimensionless Higgs field strength tp 
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in terms of which, the effective potential, Eq. (|l|), be- 
comes 



V eS ( V )=ST 
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For simplicity, we work in the approximation where the 
energy densities of both the unbroken and broken phases 
are degenerate. This happens for a value of — 2. In 
this approximation, the phase transition is described by 
a one-dimensional solution for the Higgs field, called the 
kink, which separates the two phases. This is given by 



<p(x) = 1 + tanh(ir) , 
where the dimensionless position coordinate x is 

S T 
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The parameter V2A/(<5 T) represents the width of the 
domain wall H] . It can also be checked that this param- 
eter becomes smaller in the presence of hypermagnetic 
fields. 

In terms of the kink solution we can see that x — — 00 
represents the region outside the bubble, that is the re- 
gion in the symmetric phase. Conversely, for x = +00, 
the system is inside the bubble, that is in the broken 
phase. The kink wall propagates with a velocity deter- 
mined by its interactions with the surrounding plasma. 
This velocity can be anywhere between 0.1-0.9 the speed 
of light @. 

Given the above background scalar field, the problem 
of fermion reflection and transmission through the do- 
main wall in the presence of an external magnetic field 



can be cast in terms of solving the Dirac equation. In 
the wall's rest frame, this can be written as 
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£<p{x) *(x) = 0, 



(7) 



where £ = 2m /ran, m and run being the fermion and 
Higgs masses at zero temperature, respectively and \I/ 
the fermion wave function. Also, = (0, A) is the 
four-vector potential related to the external magnetic 
field by B = V x A. In the symmetric phase, A = Ay 
corresponds to the hypermagnetic field vector potential 
whereas in the broken phase, where only the Maxwell 
projection is unscreened |^||, A = A em corresponds to 
the ordinary photon vector potential, a is the fermion 
coupling to the external field; a — ff'Y/2 in the symmet- 
ric phase with Y the fermion's hypercharge and g' the 
U(l)y coupling constant; a = e, the electric charge, in 
the broken phase. 

In the absence of a background magnetic field, the so- 
lution to Eq. (0) when considering the bubble wall as a 
planar interface has been found in Ref. |l6). Let ^(x) 
be a solution with a definite helicity (h = L, R) to Eq. (Q) 
in the absence of an external magnetic field -since we will 
be interested in fermion solutions in the symmetric phase 
where particles are massless, helicity is a good quantum 
number. It can easily be shown that in the presence of 
a magnetic field, the solution to Eq. (Q) with definite 
helicity is given by 



vMx)=^(x)e-^ Yh / XAWx ' 



(8) 



where is the hypercharge of the corresponding helicity 
mode. 

Since the presence of the hypermagnetic field shows 
up as a phase in the fermion wave function, it follows 
that its effect will dissapear when considering a process 
that involves the square of a single fermion wave function. 
However, the effect will be revealed during an interference 
process. Since this requires that at least two particles be 
present within a spatial volume of order A 3 per unit time, 
where A is the de Broglie fermion wave length, we expect 
that these processes are supressed with respect to those 
involving a single particle. We will come back later to 
estimate the probability for such processes. For the time 
being, let us consider a situation in which the fermion 
amplitude at a space point a distance I away from the 
the bubble wall receives contributions from fermions com- 
ing from the symmetric phase which are reflected from 
two points on the bubble wall separated by a distance 
d. For simplicity, we take d <^ I, in this case, we can 
consider that the fermion propagation is approximately 
directed along the direction perpendicular to the bubble 
wall. Since the coupling of the given helicity mode with 
the hypermagnetic field is chiral, the square of the ampli- 
tude representing the interference of left-handed fermions 
\M L \ 2 will differ from the square of the amplitude for 
right-handed interfering antifermions |-Mr| and conse- 
quently, an axial asymmetry A is built in front of the 
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bubble wall. Within the above approximations, it is easy 
to show that the explcit expression for the axial, CP vi- 
olating asymmetry is given by 



where 



A=\M L \ 2 - \M R \ 2 

= ^sin(T + /2)sin(T_/2), 



T± = |(Y fi ±Y i )e. 
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1Z is the fermion reflection coefficient 16 , common to left 
and right-handed modes and depends on the parameter 
£ and the fermion energy, being equal to 1 for energies 
below twice the zero temperature mass of the fermion. O 
is the hypermagnetic field flux through the area defined 
by the two trajectories of the interfering fermions. Taking 
the case of a constant hypermagnetic field perpendicular 
to this area, 6 = B Y ld/2, with By = |B Y |. Form Eq. (g) 
we see that the asymmetry disappears for a vanishing 
hypermagnetic field. 

We now follow Ref. Jljj and construct the net axial 
charge flux T in front of the expanding bubble wall in the 
symmetric phase. This flux receives contributions both 
from the reflected fermions within the symmetric phase 
and from those fermions transmitted from the bubble's 
interior. From CPT conservation, fermion transmission 
from the broken to the symmetric phase is related to 
fermion reflection in the symmetric phase. Therefore, 
the explicit expression for T is given by p7[ 
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where the bubble's wall velocity in the fluid frame is u 
and 7 = — u 2 is the Lorentz factor, pi, pt are the 

longitudinal and transverse components of the fermion's 
momentum and f s , f b are the fermion equilibrium fluxes 
(neglecting Fermi blocking factors) in the symmetric and 
broken phases, respectively 
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Since there is no net fermion number flux through the 
wall, the axial charge T and hypercharge Ty fluxes are 
trivially related by Ty = T /A. 

The rate of baryon number density production is given 
in terms of the rate of baryon number violation per unit 
volume Tb and the partial derivative of the free energy 
with respect to baryon number dF/dB by (To) 



(12) 



The quantity dF/dB corresponds to the baryon num- 
ber chemical potential ps and represents the force push- 
ing the universe towards its equilibrium baryon number 



value. Only those processes which happen fast enough 
with respect to the time that the reflected fermions in the 
symmetric phase (or those that passed from the broken 
to the symmetric phase) spend before being retaken by 
the expanding wall, will contribute to drive baryon num- 
ber towards equilibrium. This time is called the trans- 



port time t and estimates show 12 that it is of order 
t ~ 100/T. Since the only other interactions, besides 
SU(3)xU(l) gauge boson exchange and family-diagonal 
SU(2) gauge boson exchange, with a large enough rate 
are the top quark Yukawa interactions (given the large 
top Yukawa coupling), then the axial asymmetry must 
reside basically in the form of top quarks. Enforcing this 
condition together with that of havin g in itially zero net 
baryon and lepton numbers one finds |12| 
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where py is the hypercharge density, n is the number of 
scalar doublets in the theory and the minus sign comes 
from the assumption that the net axial asymmetry is in 
the form of a right-handed top number. For the SM, 
n = 1. 

It is now straightforward to integrate Eq. (|l^) to get 
the net baryon number density in terms of the axial flux, 
with the result pi] 



r B t 



(14) 



Equation (14) has to be corrected by recalling that the 
axial asymmetry, Eq. (^|), is built from two- fermion inter- 
fering processes. We estimate the probability r\ of finding 
a second fermion in the trajectory of the first one by as- 
suming that there is no initial correlation between these 
two particles. The probability to find the second one is 
thus given by the ratio of the number of particles per unit 
area and time that bounce off the bubble wall making an 
angle 9 within a solid angle subtended by a linear dimen- 
sion on the order of the particles wave length (A ~ T _1 ) 
to the incoming flux. 9 is such that 



cos 9 — 



I 



(15) 



Taking p, dp < 
r\ is given by 

V 



T and neglecting Fermi blocking factors, 



A 2 I d/2 



3ttC(3) [ Z 2 + (rf/2) 2 ] 2 1 + , 



(16) 



where C is the Riemann zeta function and the last factor 
comes from a simple Fermi-Dirac distribution evaluated 
at E = T. 

Notice that in the wall's rest frame, the interfer- 
ence point cannot be further apart than the particle's 
mean free path A m fp(~ (1 — 10)/T @), since other- 
wise fermions will rethermalize and no net baryon num- 
ber could be produced. Accounting for this probability 
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factor, the net baryon to entropy ratio produced during 
the EWPT can be written as 



Pb/s = VPb/s 



= v- 



270«a w tT 



12 



~uT 2 ' 



(17) 



where we have used that for the EWPT epoch, the en- 
tropy density is given by s — 2n 2 g+T 3 /45 with the effec- 
tive number of degrees of freedom g* ~ 107 and that the 
rate of baryon number violation per unit volume in the 
symmetric phase is given by 



T B = 3Ka w T 4 : 



(18) 



with aw = <? 2 /47r, g the SU(2) coupling constant and 

K ~ 1. 

Equation Jl7j ) is valid when the rate of baryon num- 
ber violation is much smaller than the rate at which the 
fcrmions are recaught by the expanding wall, which in 
turn sets a lower limit for the wall's velocity given by 



ina,.,TT <C u . 



(19) 



For definitiveness, we set k = 1, mu = T = 100 GeV 
and d — 1/2. (For the SU(2) and U(l)y couplings, we 
use g = 0.637 and g' = 0.344, respectively @). We find 
that for a wall velocity u — 0.1c 



p fl /s = (3-6) x HT 11 
whereas for u — 0.6c 



p B /s = (0.7-1.3) x 10 
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(20) 
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for I ~ 9/T and By = (0.3 - 0.5) T 2 . These are large 
but still not ruled out values for the background field 
strength ||, compatible with a Higgs mass of order 100 
GeV and a phase transition analog to a type I supercon- 
ductor ||. Recall that the experimental value based on 
nucleosynthesis calculations jl9| is 



p B /s = (0.1 - 1.4) x 10 



-in 



(22) 



which means that the estimates based on the present 
analysis are within the experimental limits, at least for 
slowly expanding bubbles. 

In conclusion, we have shown that in the presence of 
strong, large scale primordial hypermagnetic fields, it is 
possible to generate a large amount of CP violation that 
combined with a stronger first order EWPT -also pro- 
duced by the hypermagnetic fields- could account for the 
observed baryon number to entropy ratio within the SM. 
The fact that fermions couple chirally to background hy- 
permagnetic fields in the symmetric phase makes it pos- 
sible to build a CP violating asymmetry by considering 
two fermion interfering processes in an equivalent way to 
the Bohm-Aharanov effect. This asymmetry is converted 
into baryon number by sphaleron induced processes in 



the symmetric phase and preserved when these fermions 
are recaught by the expanding bubble wall. 

Given that some of the parameters describing the dy- 
namics of the EWPT are not very precisely determined, 
it is clear that work on that direction is necessary and 
the possibility that baryogenesis could be realized within 
the SM could certainly stimulate this kind of work. 

Support for this work has been received in part by 
CONACyT-Mexico under grant numbers 32279-E and 
29273-E. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 
[7] 



[9] 

[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 

[17] 
[18] 
[19] 



A. D. Sakharov, Pis'ma Zh. Eksp. Teor. Fiz. 5, 32 (1967) 
[JETP Lett. 5, 24 (1967)]. 

For a recent review on the subject see M. Trodden, Rev. 
Mod. Phys. 71, 1463 (1999). 

K. Kajantie, M. Laine, K. Rummukainen and M. Sha- 
posnikov, Nucl. Phys. B 466, 189 (1996). 
M. Dine, "Baryogenesis: Electroweak and otherwise", 
Proceedings of the 1994 TASI summer school CP vio- 
lation and the limits of the standard model, 507-548. 
M. Giovannini and M. E. Shaposhnikov, Phys. Rev. D 
57, 2186 (1998). 

P. Elmfors, K. Enqvist and K. Kainulainen, Phys. Lett. 
B 440, 269 (1998). 

V. Skalozub and V. Demichik, "Can baryogenesis sur- 
vive in the standard model due to strong hypermagnetic 
field?" , |hep-ph/9909550| . 

For recent reviews on the origin, evolution and some 
cosmological consequences of primordial magnetic fields 
see: K. Enqvist, Int. J. Mod. Phys. D, 331 (1998); R. 
Maartens, "Cosmological magnetic fields", International 
Confe rence on Gravitation and Cosmology, India, Jan. 
2000, astro-ph/0007352 and references therein. 



R. Durrer, P. G. Ferreira and T. Kahniashvili, Phys. Rev. 

D 61, 043001 (2000); K. Jedamzik, V. Katalinic and A. 

V. Olinto, Phys. Rev. Lett. 85, 700 (2000). 

M. Dine, O. Lechtenfield, B. Sakita, W. Fischel and J. 

Polchinski, Nucl. Phys. B 342, 381 (1990). 

A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. 

B 263, 86 (1991). 

A. E. Nelson, D. B. Kaplan and A. G. Cohen, Nucl. Phys. 
B 373, 453 (1992). 

N. Turok, Phys. Rev. Lett. 68, 1803 (1992). 

B. H. Liu, L. McLerran and N. Turok, Phys. Rev. D 46, 
2668 (1992). 

For a recent discussion on the dynamics of EWPT, see A. 
Megevand, "Development of the electroweak phase tran- 
sition and baryogenesis", hep-ph/0006177. 
A. Ayala, J. Jalilian-Marian, L. McLerran and A. P. Vis- 
cher, Phys. Rev. D 49, 5559 (1994). 
E. Torrente-Lujan, Phys. Rev. D 60, 085003 (1999). 
M. E. Carrington, Phys. Rev. D 45, 2933 (1992). 
E. Kolb and M. Turner, "The early universe" (Addison- 
Wsley, New York, 1990). 



4 



